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ABSTRACT 

The low-mass X-ray binary (LMXB) system FIRST J102347.6+003841 hosts a newly born millisec- 
ond pulsar (MSP) PSR J1023+0038 that was revealed as the first and only known rotation-powered 
MSP in a quiescent LMXB. While the system is shown to have an accretion disk before 2002, it remains 
unclear how the accretion disk has been removed in order to reveal the radio pulsation in 2007. In 
this Letter, we report the discovery of 7-rays spatially consistent with FIRST J102347. 6+003841, at a 
significance of 7 standard deviations, using data obtained by the Fermi Gamma-ray Space Telescope. 
The 7-ray spectrum can be described by a power law (PL) with a photon index of 2.9±0.2, resulting 
in an energy flux above 200 MeV of (5.5 ± 0.9) x 10 _12 ergcm -2 s _1 . The 7-rays likely originate from 
the MSP PSR J1023+0038, but also possibly from an intrabinary shock between the pulsar and its 
companion star. To complement the 7-ray study, we also re-investigate the XMM-Newton data taken 
in 2004 and 2008. Our X-ray spectral analysis suggests that a broken PL with two distinct photon 
indices describes the X-ray data significantly better than a single PL. This indicates that there exists 
two components and that both components appear to vary with the orbital phase. The evidence for 7- 
ray emission conforms with a recent suggestion that 7-rays from PSR J1023+0038 may be responsible 
for ejecting the disk material out of the system. 

Subject headings: gamma rays: stars — Pulsars: individual (PSR 1023+0038) — X-rays: binaries 



1. INTRODUCTION 

Radio millisecond pulsars (MSPs) are rotating neutron 
stars (NS) that have been spun up via the transfer of an- 
gular momentum throu gh accretion in low-mass X-ray bi- 
nary (LMXB) systems (Bhattac harva fc van den Heuvell 
1991). Recently, a total of 9 7-ray emitting 
MSPs have been identified thr ough their 7-ray pulsa- 
tions (lAbdo et alJI200^[2010ah . 

Theoretical models have long suggested that an 
accretion-powered MSP in a LMXB will turn on as 
a rota tion-powered MSP when the system is in quies- 
cence (|Alpar et alj 119821) . How exactly the transition 
happens is not completely understood, however, it is 
widely believed that radio MSPs can only turn on af- 
ter the accretion disk is removed. Sugg ested mech- 
anism s include the pulsar wind ablation (jWang et al.1 
2009), heating associated with the d eposition of 7-rays 
from t he MSP dTakata et a!ll2010aD . and the propeller 
effect (jRomanova et al.M2009D ■ 

PSR J1023+0038 is the first and only known 
rotation-powered MSP in a quiescent LMXB, 
namely FIRST J102347.67+003841.2 (hereafter 
J1023) J1023 was identified as a LMXB in 
2006 (jHomer et al.1 I2006D and the radio M SP was 
found subsequently (Archibald e t al.l 12009ft 
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2002 dWang et all I2009D and th e disk has since 
disappeared ([Archibald et al.1 I2009D ; rad io pulsation 
was found in 2007 (jArchibald et all 120091) . Therefore 
PSR J1023+0038 is considered as a newly born MSP, 
representing the long sought-after missing link of a 
rotation-powered MSP descended from a LMXB. 
The spin dow n power of PSR J10 23+0038 (L sd < 



3 x 10 35 erg s^ 1 : lArchibald et al.ir2009D is relatively high 



compared to other 7-ray MSPs (|Abdo et al.ll2009afl . A 
small fraction of this power would suffice to generate de- 
tectable 7-rays from PSR J1023+0038. Given the recent 
evidence of an accretion disk before 2002 that later dis- 
appeared, it is of great importance to probe the energy 
source that facilitated the dissolutio n of the disk. The 
detection of X-rays from the system (Ho mer et al.ll200"6T) 
provides a hint that high-energy processes are ongoing, 
which is further strengthened by the reported X-ray pul- 
sation f rom PSR J1023+0038 and orbital variability in 
X-rays (|Archibald et al.ll2010D . These results prompted 
us to search for 7-rays from the system and to study its 
X-ray properties in more detail. 

2. GAMMA-RAY OBSERVATIONS AND ANALYSIS 
RESULTS 

The Large Area Telescope (LAT) aboard the 
Fermi Gamma-ray Space Telescope is able to de- 
tect 7-ray s with energies b etween ~20 MeV and 
>300 GeV (jAtwood et al.ll2009D . Data used in this work 
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were obtained between 2008 August 4 and 2010 July 14 
that are available at the Fermi Science Support CenteiQ. 
We used the Fermi Science Tools v9rl5p2 package to re- 
duce and analyze the data in the vicinity of J1023. Only 
those data that passed the most stringent photon selec- 
tion criteria (i.e. the "diffuse" class) were used. To re- 
duce the contamination from Earth albedo 7-rays, we 
excluded events with zenith angles greater than 105°. 
The instrument response functions "P6_V3_DIFFUSE" 
recommended for analysis of the "diffuse" class events 
were used. We chose 200 MeV as the lower energy cut 
to include sufficient source photons at low energies while 
reducing the contamination of background photons that 
dominates at low energies. Therefore, we used events 
with energies between 200 MeV and 20 GeV in the like- 
lihood analysis. 

We carried out an unbinned maximum-likelihood spec- 
tral analysis (gtlike) of the circular region of 15° radius 
centered on the 7-ray position (see below). We sub- 
tracted the background contribution by including the 
Galactic diffuse model (gll_iem_v02.fit) and the isotropic 
background (isotropic_iem_v02.txt), as well a s all sources 
in the first Fermi/LAT catalog (1FGL; lAbdo et all 
l2010bt l within the circular region of 25° radius around 
the 7-ray position. We assumed a power law (PL) spec- 
trum for all the 52 1FGL sources considered. The spec- 
tral parameter values were set free for sources within 10° 
from J1023. 

Th e maximized test-statistic (TS) value (jMattox et al.l 
1996) we obtained for the pulsar position is 50, corre- 
sponding to a detection significance of 7a. The position 
of the 7-ray source is estimated by gtfindsrc to be at 
R.A. (J2000) = 155?92 and Dec. (J2000) = 0?72 with 
statistical uncertainty 0?08 (0?2) at the 68%(95%) con- 
fidence level, which is consistent with the position of 
J 1023. The systematic u ncertainty is estimated to be 
<40% (jAbdo et al.ll2009bfi . We used gttsmap to obtain 
the TS map of the 5°x5° region centered on the best-fit 
7-ray position, as shown in Fig. [TJ 

To investigate why this 7-ray source was not present 
in the 1FGL catalog, we divided the whole dataset into 
two: the first and second year data, respectively, and 
performed the likelihood analysis for each of them sepa- 
rately. Both TS values drop to 22-24, which are about 
half of the TS value derived using the whole dataset 
and below the threshold of TS=25 to be included in the 
1FGL catalog (jAbdo et al.l l2010bf l. This suggests that 
the source we found is significant (i.e., >5<t) with two- 
year data. The spectral parameters from each dataset 
are also consistent with the values reported below. 

We then fit the 7-ray spectrum with a single PL, result- 
ing in a photon index of 2. 9 ±0.2 and an integrated energy 
flux above 200 MeV of (5.5 ± 0.9) x 10~ 12 erg cm~ 2 s -1 . 
We divided the 200 MeV-20 GeV 7-rays into five en- 
ergy bins of logarithmically equal bandwidths and recon- 
structed the flux using gtlike for each band separately, 
assuming a PL model with T 7 =2.9 within each bin. No 
7-ray source was needed at the J1023 position in the 
three bins above 1.3 GeV (derived TS values <5) in the 
likelihood analysis, indicating a cut-off at ~1 GeV. We 
therefore attempted to fit the 200 MeV-20 GeV spectrum 
with a PL with an exponential cut-off (PLE). We found 

1 http: //fermi .gsfc .nasa.gov/ssc/ 



that PLE also well describes the spectrum, i.e. both PL 
and PLE models gave the same TS value of 50. However, 
it should be noted that a spectral cutoff is statistically 
not required. Due to the low photon statistics (~280 
modeled photons from J1023), the spectral parameters 
of the PLE model cannot be well constrained simultane- 
ously. Motivated by the possible magnetospheric origin 
of the 7-rays, T-y and E c are estimated by fixing the other 
one (at its mean value) to be 1.9±0.3 and 700±220 MeV. 
We fixed L 7 at 1.5 to 2.4 (with steps of 0.1) while letting 
the normalization and E c free, and found that E c was 
only well constrained when T 7 = 1.7 to 2.1, consistent 
with the above result. 

At the distance of 1.3 kpc, the 7-ray luminosity 
(above 200 MeV) is (1.1 ± 0.2) x 10 33 erg s" 1 . As- 
suming that the 7-rays come from PSR J1023+0038 
(see Sect. 131, th e pulsar spin-down l uminosity E < 
3 x 10 35 ergs s" 1 ([Archibald et al.ll2009D implies a 7-ray 
conversion efficiency of only >0.3%. Such 7-ray luminos- 
ity and 7-ray efficiency are among the smallest of 7-ray 
MSPs. Moreover, assuming the PLE model is robust, 
the cutoff energy of ^700 MeV is the lowest among all 
7-ray MSPs. The spectral properties of the 7-ray emis- 
sion from J1023 are summarized in Table [TJ 

A search for pulsation of the >10 MeV 7-rays within 
a 1° radius region around the 7-ray position was per- 
formed. We did not find any significant pulsed detection 
at or close to the spin-period of the pulsar nor any indica- 
tion of 7-ray variability related to the orbital modulation. 
Even so, given the low photon statistics (of just over 700 
photons), this result does not preclude any 7-ray pulsa- 
tion. 

We also performed a long-term temporal analysis of 
J 1023, in which the 2-year data were binned in 3-month 
periods. No significant 7-ray variability was found, indi- 
cating that the object is stable (down to 3-month peri- 
ods) in radiating 7-rays. 

3. A MODEL OF 7-RAY EMISSION FROM PSR J1023+0038 

The 7-rays from J1023 originate either from the pul- 
sar magnetosphere or a shock where material over- 
flowing from the companion interacts with the pul- 
sar wind. In the pulsar wind scenario, 7-rays may 
be generated as synchrotron radiation of the electrons 
and p ositrons accelerated in the shock (jTakata fc Taaml 
2009). An important theoretical uncertainty is how 
a particle kinetic-dominated flow can be formed near 
the pulsar ([Kirk fc Ski r] [200l lAronsl[2008l) . In 

addition, it is also difficult to explain the steep PL 
spectrum (r 7 ^3) in the synchrotron model. On the 
other hand, if the pulsar magnetosphere is sufficiently 
clear of matt er, GeV 7-ray photo ns can be produced 
in a slot- gap ([Harding et al.l 20051 ) or an outer-gap ac- 
celerator ( Cheng et al.l H986aHbh . The curvature radia- 
tion from a gap accelerator is expected to have a PLE 
shape. The reported X -ray pulses from PSR J1023+0038 
([Archibald et al.ll2010f) suggests that the observed 7-rays 
are also produced in the pulsar magnetosphere. More 
accumulated source photons in the future should help to 
answer whether there exists 7-ray pulsation or variability 
on the time scale of the orbital period, in turn helping to 
distinguish these two scenarios for the observed 7-rays. 

Applying the outer gap model, the 7-ray luminosity is 
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approximately described by 

£7 ~ f 3 L s d, (1) 

where / is the fractional gap thickness, i.e. the ratio of 
the gap thickness to the radius of the light cylinder and 
L s d = 4:(2n) 4 B 2 R 6 /6c 3 P 4 the pulsar spin down power 
with B corresponding to the surface magnetic field, R 
the stellar radius, and P the rotation period. The frac- 
tional gap thickness is determined by the pair-creation 
condition between 7-rays emitted in the outer gap and 
X-rays from the NS surface, i.e., E 1 E S: x — (m-eC 2 ) 2 , 
where E 1 ~ 1 GeV and E Si x are the typical energy 
of the 7-rays and the X-rays, respectively. Assuming 
E s . x = 30 eV ( s ee Sect . H and Table [2]), and follow- 
ing rTakata et all (|2010bf ) we find f - 0.5. The 7-ray 
luminosity is estimated as L 1 ~ 4 x 10 34 erg s _1 using 
^sd = 3 x 10 35 erg s _1 , consistent with the required lumi- 
nosity from an irradia ting source to explain the heatin g 
of the companion star ([Thorstensen fc ArmstronsJl2005l ). 

We may estimate the 7-ray flux measured on the 
Earth as F 1 ~ L 1 /8Q,d 2 ~ 10 -9 erg/cm 2 s, where L 7 = 
3 x 10 34 erg s _1 , the solid angle 5fl = 2, and the dis- 
tance d = 1.3 kpc are assumed. This flux is higher than 
the observed F 1 ~ 5 x 10 -12 erg cm -2 s _1 . We note, 
however, that the 7-ray flux depends on the viewing ge- 
ometry, because the intensity varies over the 7-ray beam. 
We carried out a simple 3-D calculati on using curvature 
radiation in the ou ter gap model (c.f. iWang et al.ll2010t 
iTakata et all l2007[ ) and found that F 7 increases from 
~5xl0 -13 erg cm~ 2 s _1 to ^10~ n erg cm~ 2 s _1 when 
the viewing angle, i.e. the angle between the rotational 
axis of the pulsar and the line of sight, £ , increases from 
-34° to -53° ( see I Archibald et al.ll2009l) . well consistent 
with the observed 7-ray flux. 

4. RE-ANALYSIS OF THE X-RAY DATA 

XMM-Newton observations of J1023 were conveyed 
on 12 May 2004 (hereafter we refer this observation as 
XMM1) with all the EPIC cameras operated in full-frame 
mode, and on 26 November 2008 (hereafter XMM2) with 
the MOS 1/2 CCDs operated in full-frame mode and the 
PN camera operated in timing mode with a resolution of 
0.03 ms. 

Archibald et al. (2010) reported an X-ray analy- 
sis of J1023, finding a possible X-ray pulsation from 
PSR J1023+0038 and modulation on the orbital period 
of the binary. They interpreted the latter as coming from 
an intrabinary shock, and suggested a composite spectral 
model with a PL plus a possible thermal component. Al- 
though this model results in an acceptable goodness-of- 
fit, we found that their analysis is apparently incomplete. 
The X-ray spectrum from the pulsar magnetosphere is 
typically different from the shock emission (cf. Hui & 
Becker 2006, 2007, 2008), and the combined contribu- 
tion from the two emission regions cannot be generally 
described by a single PL. To complement their results, we 
report our independent analysis of these XMM-Newton 
observations. 

As PN data obtained from XMM2 are collapsed into 
a one-dimensional row, they are not suitable for spec- 
troscopy. After filtering for the high sky background and 
events affected by bad pixels, the effective exposures are 
14.9 (MOS 1/2) and 11.6 kiloseconds (PN) for XMM1, 



and 33.5 kiloseconds (MOS1/2) for XMM2. 

To avoid contamination from a nearby X-ray source, we 
extracted the energy spectrum of J1023 from the circles 
with radius 25" centered on the radio timing position 
in all datasets, which corresponds to an encircled energy 
fraction ~ 80%. The background spectra were sampled 
from nearby low-count circular regions of radius 40" in 
the corresponding cameras. We grouped each spectrum 
dynamically to obtain the same signal-to-noise ratio in 
each dataset. 

First, we fit the XMM1 and XMM2 data separately 
with various single-component models. Only the PL 
model provides a good description of the data in both 
XMM10 and XMM2 (see Table [2]). With no indication 
of spectral and flux variability between data taken in 
XMM1 and XMM2, we combined both datasets for a 
constraining spectral analysis. 

We found that the single PL model provides an ac- 
ceptable description of the combined data with the spec- 
tral parameters as in Archibald et al. (2010). However, 
we found systematic deviations in this model, indicating 
that additional componcnt(s) might be required. Mo- 
tivated by the possible presence of both pulsed emis- 
sion and emission dependent on the orbital phase of 
the binary, we fit the spectrum with a broken power 
law (BKPL). The goodness-of-fit was found to be im- 
proved significantly (x 2 = 0.75 for 75 d.o.f.) and no 
systematic fitting residuals were found (see Fig. [2]), in 
contrast with the best- fit model (i.e. s i ngle P L model) 
presented in Fig. 1 of lArchibald et all (|20100 . To sta- 
tistically address the improvement for the spectral fits 
from the single PL to the broken PL, we used an F- 
test and found that the p- value is 7.7xl0~ 6 . There- 
fore, the additional parameters are required statistically. 
This behavior was found in both epochs (XMM1 as 
well as XMM2). The break energy was found to be 
E b = 1.84±°,; 2 g keV. The spectrum was found to be 
steeper (i.e. T x = 1.75±o.h) at E < E b than that in 
the hard band (i.e. T 2 X = 1.07±° °£). The unabsorbed 
flux was found to be fx = 6.3^J'g x 10 -13 ergs cm~ 2 s _1 
(0.01-10 keV). At a distance of d ~ 1 kpc, the isotropic 
X-ray luminosity is L x = 7.5t 4 9 x 10 31 ergs s _1 (0.01- 
10 keV), corresponding to an X-ray conversion efficiency 
of > 0.03% of the spin-down luminosity. 

As thermal X-rays from the NS surface may also con- 
tribute, we estimated the thermal contribution besides 
the NMHA model with the mass and radius of the NS 
fixed at Mns = 1-4 M© and i?NS = 10 km. Neverthe- 
less, this additional thermal model was only required at 
a confidence level of < 50%. Furthermore, with nor- 
malization as a free parameter, an unreasonably small 
source distance of ~ 30 pc was inferred. To constrain 
the possible thermal contribution from the stellar sur- 
face, we fixed the normalization corresponding to 1 kpc. 
This results in the la upper limits of the temperature 
and the thermal flux contribution to be kT < 32.76 eV 
and fx < 7.3 x 10~ 14 ergs cm" 2 s" 1 (0.01-10 keV), re- 

2 Circles of radius 35" were used infArchibald et al. (2010), that 
probably contain a larger amount of contamination from the PSF 
wing of that nearby X-ray source. 

3 Our results for XMM1 are consistent with those reported by 
Homer et al. (2006). 
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spectively. 

The fact that more than one PL index is required to 
explain the observed X-ray spectrum strongly suggests 
an additional X-ray contribution from this system be- 
sides the magnetospheric radiation. As possible orbital 
modulation in X-rays has been identified, the additional 
non-th ermal component may come from the intrabinary 
shock ((Archibald et al.lfeOlOl) . Furthermore, our spec- 
tral analysis suggests that the emission below and above 
the break energy (—2 ke V) might have different origins. 
lArchibald et alJ (|2010t ) reported an X-ray pulsation for 
photons at 0.25-2.5 keV. Whether the emission above 
—2 keV is pulsed remains unclear. 

We performed an energy-resolved analysis of the X- 
ray flux as a function of orbital phase. According to 
the resultant break energy, we divided the energy range 
(0.3-10 keV) suitable for timing analysis into two bands: 
0.3-2 keV (soft) and 2-10 keV (hard). We found no sig- 
nificant variability at the orbital period in both band s 
from XMM1 data, consistent with iHomer et al.l (|2006l ). 
This is likely due to low photon statistics, as well as insuf- 
ficient orbital coverage of -91% for MOS1/2 and -82% 
for PN data. We therefore focused on the XMM2 data in 
this study. The orbital light curves of the soft and hard 
bands are shown in Fig. [31 Both soft and hard X-rays 
contribute to the orbital modulation and the troughs and 
peaks in the orbital light curves in these two bands occur 
at the same orbital phases (around 0.2 and 0.6 respec- 
tively). Using a x 2 -test, the significance for a flux modu- 
lation over the observed orbit in the soft and hard bands 
was found to be —99.1% and —85.5%, respectively. Al- 
though the significance derived from this test is higher for 
soft X-rays than for hard X-rays, we could not determine 
the energy band for which the modulation is greater. 



5. CONCLUDING REMARKS 

We have found strong evidence for 7-ray emission from 
a LMXB system, FIRST J102347.67+003841.2. The 
7-rays may originate from PSR J1023+0038, or possi- 
bly from the intrabinary shock between the pulsar and 
its companion star. Given the observed steep spec- 
trum (r 7 —3), synchrotron emission from the shock 
is not favored. We have applied an out er gap model 
to expl ain the observed 7-rays. Recently, iTakata et al.l 
(|2010af) suggest that such 7-rays from a newly-born MSP 
may be responsible for the disk clearance in such sys- 
tems. Our discovery of 7-rays, once confirmed, from 
FIRST J102347.67+003841.2 conform with their predic- 
tions. 

The spectral analysis of the XMM data taken in 2004 
and 2008 suggests that a broken PL well describes the 
X-ray data, with the lowest systematics, compared to a 
single PL and other composite models tested. Our X-ray 
spectral analysis suggests that there exists two compo- 
nents and that both components appear to vary with 
th e orbital phase, suppor ting the X-ray variability study 
bv l Archibald et af] (|2010D . that suggested more than one 
mechanism contributes to the observed X-ray flux from 



J1023. 

In such a system that recently transited from the 
LMXB phase to the radio MSP phase, we speculate that 
the current state may not be permanent as an accre- 
tion disk may reform about the NS in PSR J1023+0038 
switching off the radio MSP. A cycle of on and off states 
of pulsar activity could repeat in the near future. Contin- 
uous monitoring in all wavelengths of this unique source 
will test this idea and help us to comprehend the on-going 
evolution of FIRST J102347.67+003841.2. 

We acknowledge the use of data and software facil- 
ities from the FSSC, managed by the HE AS ARC at 
the Goddard Space Flight Center. CYH is supported 
by research fund of Chungnam National University in 
2010. AKHK and LCCL are supported partly by the Na- 
tional Science Council of the Republic of China (Taiwan) 
through grant NSC99-2112-M-007-004-MY3 and NSC99- 
2811-M-008-057, respectively. YJY has received fund- 
ing from the European Community's Seventh Frame- 
work Programme (FP7/2007-2013) under grant agree- 
ment number ITN 215212 "Black Hole Universe". KSC 
is supported by a GRF grant of Hong Kong Government 
under HKU700908P. 



REFERENCES 



Abdo, A. A., et al. (Fermi/LAT Collaboration) 2009a, Science, 
325, 848 

Abdo, A. A., et al. (Fermi/LAT Collaboration) 2009b, ApJS, 183, 
46 

Abdo, A. A., et al. (Fermi/LAT Collaboration) 2010a, ApJ, 712, 
957 

Abdo, A. A., et al. (Fermi/LAT Collaboration) 2010b, ApJS, 188, 
405 

Alpar, M. A., Cheng, A. F.. Ruderman, M. A., & Shaham, J. 

1982, Nature, 300, 728 
Archibald, A. M., et al. 2009, Science, 324, 1411 
Archibald, A. M., et al. 2010, ApJ, 722, 88 
Arons, J. 2008, AIP Conf. Proa, 983, 200 

Atwood, W. B., et al. (Fermi/LAT Collaboration) 2009, ApJ, 
697, 1071 

Bhattacharya, D., van den Heuvel, E. P. J. 1991, Phys. Rep., 203, 
1 

Cheng, K. S., Ho, C, & Ruderman, M. 1986a, ApJ, 300, 500 
Cheng, K. S., Ho, C, & Ruderman, M. 1986b, ApJ, 300, 522 
Harding, A. K., Usov, V. V., & Muslimov, A. G. 2005, ApJ, 622, 
531 

Homer et al. 2006, AJ, 131, 562 

Hui, C. Y. & Becker, W. 2006, A&A, 448, L13 

Hui, C. Y. & Becker, W. 2007, A&A, 470, 965 

Hui, C. Y. & Becker, W. 2008, A&A, 486, 485 

Kirk, J. G. & Skjaeraasen, O. 2003, ApJ, 591, 366 

Mattox, J. R., et al. 1996, ApJ, 461, 396 2002, A&A, 387, 993 

Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & 

Lovelace, R. V. E. 2009, MNRAS, 399, 1802 
Takata, J. & Taam, R. E. 2009, ApJ, 702, 100 
Takata, J., Chang, H.-K., & Cheng, K. S. 2007, ApJ, 656, 1044 
Takata, J., Cheng, K. S., & Taam, R. E. 2010a, ApJ, 723, L68 
Takata, J., Wang, Y., & Cheng, K. S. 2010b, ApJ, 715, 1318 
Thorstensen, J. R., & Armstrong, E. 2005, AJ, 130, 759 
Wang, Y., Takata, J., & Cheng, K. S. 2010, ApJ, 720, 178 
Wang, Z., Archibald, A. M., Thorstensen, J. R., Kaspi, V. M., 

Lorimer, D. R., Stairs, I., & Ransom, S. 2009, ApJ, 703, 2017 



5 



Table 1 

Spectral parameter values of 7-ray emission from J1023. 



Model 3 


Photon flux b (>200 MeV) 


Energy flux (>200 MeV) 


Photon Index 


Cutoff energy 


7-ray luminosity (>200 MeV) 




(10- 9 cm- 2 s- 1 ) 


(ltr 12 erg cm- 2 s" 1 ) 


r 7 


(MeV) 


(10 33 erg s- 1 ) 


PL 


8.2 ±1.6 


5.5 ±0.9 


2.9 ± 0.2 




1.11 ± 0.18 


PLE 


8.1 ± 1.5 


5.2 ±0.8 


1.9 d 


700 ± 230 


1.06 ±0.17 


PLE 


7.8 ±1.5 


5.1 ±0.8 


1.9 ± 0.3 


700 


1.03 ±0.17 



a PL=power law model; PLE=power law with an exponential cut-off model 

b All the quoted errors are statistical and la for one parameter of interest. 

c The pulsar distance is taken as 1.3 kpc. 

d Model parameters without quoted errors are fixed at the value given. 



Table 2 

Spectral parameter values of X-ray emission from J1023. 



Model 3 


Epoch 


n H b 
(10 2() cm- 2 ) 


Photon Index c 


E b /kT 
(keV) 


Unabsorbed flux d 
(10- 13 crg cm- 2 s- 1 ) 


Xl (d.o.f.) 


PL 
PL 


12 May 2004 (XMM1) 
26 Nov 2008 (XMM2) 


<1.2 
<0.4 


1 30+ -" 4 
1.29 ± 0.04 


.../... 


4 6 +0 - 5 

4 Q+0.4 
^• M -0.3 


0.94 (40) 
1.17 (35) 


PL 

BKPL 
BKPL 
& NMHA 


XMM1 & XMM2 
XMM1 & XMM2 
XMM1 & XMM2 


<0.6 

4.4 ± 1.6 
4 7 +3.3 


1.29 ±0.03 
1.75t°"i; (E < E b )/1.07±° f 6 (E > E b ) 
1.76t°;i!(£ < E b )/im+_l X g(E > E b ) 


1.84+°^/... 
l-84t°;^/<0.033 


4 7 +0.3 
' -0.2 
6 3+ 1 - 6 

6.4±J"| (BKPL) 
< 0.73 (NMHA) 


1.00 (77) 
0.75 (75) 
0.76 (74) 



a PL— power law model; BKPL— broken power law model; NMHA— non-magnetic hydrogen atmospheric model 

b column absorption 

c All quoted errors are ler for one parameter of interest. 

d in the 0.01-10 keV range 
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Figure 1. Test-statistic (TS) map of a region of 5° x5° centered at the position of J1023 (labeled by the cross). Gamma-rays with energies 
between 200 MeV and 20 GeV were used. This map is created by moving a putative point source through a grid of locations on the sky 
and maximizing — log(likelihood) at each grid point, while stronger and well-identified sources outside the sky map are included in each fit. 
The Fermi source 1FGL J1022. 8-0115 is clearly seen in the map. The 95% con fidence- level error c ircles of the best-fit position of the 7-ray 
emission are also shown. The error circle of lFGL J1022. 8-0115 is taken from Abdo ct al. ( 2010bT). 




Figure 2. The X-ray spectra of J1023 as obtained from the observations taken at 12 May 2004 and 26 November 2008 with XMM-Ncwton 
and simultaneously fitted to an absorbed broken power-law model (upper panel) and contribution to the \ 2 fit statistic (lower panel). 
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Figure 3. The energy-resolved X-ray light curves of J1023 in the soft (0.3—2.0 keV) and hard (2.0-10.0 keV) bands, as obtained from the 
observation taken at 26 November 2008 with XMM-Ncwton. 



